Spin Splitting in Symmetrical SiGe Quantum Wells 
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Spin splitting of conduction electron states has been analyzed for all possible point symmetries of 
SiGe quantum well structures. A particular attention is paid to removal of spin degeneracy caused by 
the rotoinversion asymmetry of a (001) heterointerface between two diamond-lattice materials. The 
asymmetry is shown to result in spin splitting of both Rashba and Dresselhaus types in symmetrical 
SiGe quantum wells. Consequences of the spin splitting on spin relaxation are discussed. 



Spin properties attract the great attention in recent 
years due to attempts to realize an electronic device 
based on the spin of carriers. Conduction electrons 
are obvious candidates for such devices, particularly 
in nanostructures where electron energy spectrum and 
shape of the envelope functions can be effectively engi- 
neered by the growth design, application of electric or 
magnetic fields as well as by illumination with light. 

Various semiconductor materials are being involved 
in the spintronics activities. SiGe quantum well (QW) 
structures are among them. 1-2-3 - 4 Although bulk Si and 
Ge have an inversion center, QW structures grown from 
these materials can lack such a center and allow the spin 
splitting of the electronic subbands. 

The quantum engineering of spintronic devices is usu- 
ally focused on the Rashba spin-dependent term to the 
electron effective Hamiltonian in heterostructures. This 
contribution appears due to asymmetry of the heteropo- 
tential (the so-called Structure Inversion Asymmetry, or 
SIA) and has no relation to the properties of a bulk semi- 
conductor. 

In III-V heterostructures, there exists another spin- 
dependent contribution called the Dresselhaus term that 
appears due to Bulk Inversion Asymmetry (BIA). It is 
commonly believed that the Dresselhaus contribution is 
absent in structures grown from centrosymmetric mate- 
rials. 

In the present work we show that the Dresselhaus-like 
spin splitting is possible in heterostructures made of Si 
and Ge due to the anisotropy of chemical bonds at inter- 
faces. 

Symmetry of a (OOl)-grown interface between 
Sii-zGea; and Si can be C2 V or C±v on averaged The 
former point group describes the symmetry of an ideal 
heterointerface with the interfacial chemical bonds ly- 
ing in the same plane. A nonidcal interface containing 
monoatomic fluctuations has two kinds of flat areas with 
interfacial planes shifted with respect to each other by 
a quarter of the lattice constant. The local symmetry 
of each area is C2 V as well. However if the both kinds 
are equally distributed, the interface overall symmetry 
increases up to C^,. It follows then that the symmetry 
of a Sii-^Ge^/Si QW structure containing two interfaces 
is described by one of five point groups: i?2d or Z?2h in 
case of two ideal interfaces with odd or even number of 
monolayers between them; C2 V for a pair of ideal and 



rough interfaces; Ci V or for two non- ideal interfaces 
of the overall symmetry C^ v each, see Fig. 1. 

Two of the above-mentioned groups, and £>2h, 

contain the space inversion operation and forbid the spin 
splitting of electronic states. Three remaining groups al- 
low the spin-dependent linear- in-fc term Tlk = "lafiUakp. 
Here a a are the Pauli matrices and fc = (k x ,k v ) is the 
in-plane electron wave vector. 

In Sii-zGe^/Si QWs grown along direction z || [001] 
with low enough content of Ge in the alloy layer, the low- 
est conduction band is located near the X point of the 
Brillouin zone. Note that in the following we consider the 
electronic states attached to the X z valley because, due 
to the quantum confinement effect, its bottom lies lower 
than X x and X y valleys^ Hereafter the subscripts x, y, z 
indicate three X valleys as well as the Cartesian coordi- 
nate system with x || [100], y || [010] being the in-plane 
axes. The symmetry analysis of the electron quantum 
confined states is based on the fact that at the X point 
the bulk Bloch functions form projective representations 
of the point group D^. All five above-mentioned groups 
are subgroups of . In the group Z?2d , the Dresselhaus- 
like term 



h D (k) 



Oyky 



is an only invariant which can be constructed from the 
products a a kp. On the contrary, in the group C^ Vl the 
only invariant combination is the Rashba term 



h R (k) = a x k y 



Oyk X 



The analysis shows that both combinations ho{k) and 
h R (k) are invariants of the group C2 V , i-e. both 
Dresselhaus- and Rashba-like spin-dependent terms are 
allowed. 

In SiGe/Si QWs with high Ge content, the lowest con- 
duction band may be located at the L point. The bulk 
Bloch functions form projective representations of the 
Z?3d point group which are p-equivalent to the usual rep- 
resentations of the same group. In this case the associ- 
ated coordinate system x',y',z' is connected with the 
principal valley axis z' || [111] and the in-plane axes 
x' || [110] (perpendicular to one of the mirror-reflection 
planes) and y' || [112]. In (OOl)-grown QWs of the sym- 
metry £>2h or Z)4h, intersection of these point groups with 
D 3 h is C*2h with the inversion, so that the spin degener- 
acy is retained. On the other hand, the intersection of 
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FIG. 1: Different interface profiles and QW point symmetries. The growth direction is [001], n is a number of monoatomic 
layers. 



the rest three groups with Z? 3 d is C s . As a result, the 
effective Hamiltonian for the L-valley electrons in QWs 
of the symmetry Z?2d , or C^ v contains three linearly 
independent combinations o~ y > k x > , ay k x > and ay k y > re- 
sponsible for the spin splitting. 

SiGe/Si heterostructures grown in the (11 Indirection 
have cither or C*3 V point symmetry depending on the 
parity of monolayer numbers. The group has the in- 
version center and retains the spin degeneracy whereas 
the C3V group allows the spin splitting. The lowest con- 
duction band in (lll)-SiGe/Si QWs is located at the L 
point of the Brillouin zone. Since Cs v is a subgroup of 
Z?3d, the relevant symmetry is Cs v , and the invariant 
combination of the products o~ a kp is the Rashba term 

(7x'ky' (Jy'kx' • 

Microscopically, a linear-in-fc correction to the 
conduction-band effective Hamiltonian is given by the 
second order of the perturbation theory 



AH = 



H soHjp + HkpHso 

Er — En, 



(1) 



Here H kp and -ffso are interband k • p and spin-orbit in- 
teraction Hamiltonians, respectively. We take into con- 
sideration only the coupling of the conduction X\- and 
the valence ^-states. The symmetry properties of the 
Bloch functions at the X point of a diamond-lattice semi- 
conductor crystal coincide with those of the following 
functions 7 



Xi-states 



S = cos (2irz/a) 
Z = sin (2ixzj a) 



Xi-states 



X = sin (2irx/a) cos (27ry/a) 
Y = cos (2irx/a) sin (2iry/a) 



(3) 



where a is the lattice constant. For the bulk states X\ , X4 
in the bases (J2J, ©, the interband matrix elements can 
be presented as 



H 



kp 



k 



H 



so 



■ U 



(4) 



Here V = (h/m ) (S\p x \X), U = (S\U X \X), p is 
the momentum operator, and the pseudovector U = 



(h/4mQC 2 )'VW x p enters into the spin-orbit Hamilto- 
nian Hso — • U (W is the microscopic potential). 
The substitution of (0J into results in 



AH oc h D {k) 



1 

-1 



where the second multiplier is a 2x2 matrix related to 
the basis (0). The matrix AH. does not lead to lifting the 
degeneracy of the Xi-states in the bulk centrosymmet- 
ric material, in accordance with the general symmetry 
consideration. 

The splitting can be achieved if one takes into account 
the anisotropy of chemical bonds at the interfaces. It 
results in (5-functional contributions to H kp and -ffso of 
the form 

AH kp = Vk P 8(z - Zif), AH so = V so S(z - z if ), (5) 

where Zif is the interface coordinate, and the matrices 
Vfc P ,Vso have few linearly independent components. In 
the case of the lowest symmetry under study, C*2 V , each 
matrix is determined by four linearly-independent pa- 
rameters 



V k% 



Vso 



Plk X +P 2 ky P 1 ky+P 2 k x 

P 3 k y + P A k x P 3 k x + P 4 k y 



U\o x — Ui<Jy —UiOy + U20- x 



(6) 



(7) 



(n 



In the bases (0), ©, the parameters P n and U n 
1-^4) are purely imaginary. 

By using Eqs. and 10-0, one can show that, for 
a single interface of the C2 V -symmetry, the correction 
(2) to the conduction-band Hamiltonian linear in k x t y and 
responsible for the removal of spin degeneracy is given 
by 



AH 



c 2v 



H if S(z-Zif) 



(8) 



Here Hif is a linear combination of the spin Pauli matri- 
ces that, in the first-order approximation in the pertur- 
bations JSJ, has the form 



Hif - — 



E 
U 
Eo 



(U 3 - Ux) h D {k) + (Ui - U 2 ) h R {k) 
{Pi + P3) h D (k) + (P 2 + P^ h R {k) 



(9) 



3 



with Eg being the band gap between X\ and X4 states. 

The electron effective Hamiltonian in an ideal QW con- 
tains a sum of two contributions ® related to the left- (I) 
and right-hand-side (r) interfaces. If the QW contains 
an even number of monoatomic layers (D 2 \ 1 symmetry), 
then the corresponding parameters U l n ,U^ or P^,P£ are 
interconnected by 

[Ui, U 2 , U 3 , Ui] r = [U 3 , U 4 , Ux, U 2 ]i, 



[Pl,P 2 , Pi, Pi]r = [-P3, -Pi, -Pl,-P 



2\l: 



and if the QW contains an odd number of monoatomic 
layers (D 2< x symmetry), then 

[Ux, U 2 , Us, Ui] r = [Ux,-U 2 , Us, -Ui] u 



[Px,P2,P 3 ,Pi]r = [Pl,-P2, Ps, -Pi]l- 

From here it readily follows that, in the former case, the 
electronic states in SiGe QWs are spin-degenerate (Hi f = 
0) and, in the latter, the spin degeneracy is removed by a 
Dresselhaus-like linear-in-fc terms. The similar analysis 
can be carried out for non-ideal SiGe QW structures of 
the C 2v and C 4v symmetry. 

Spin splitting of L-electrons. Without spin, the basis 
functions for the L-point may be chosen as s-like for the 
L\ conduction-band state (the Bloch function S'), and 
as p-like (X',Y') for the Ls< valence-band states. With 
spin, they are multiplied by the spin functions \ and \. 

In a bulk diamond-lattice semiconductor, the spin- 
orbit interaction results in a splitting by some value A 
of the Ls> valence band into the (L^,L^) and L§ sub- 
bands. In the three-fold degenerate T-point, the second- 
order perturbation theory with allowance for A^O gives 
rise to Dresselhaus-like conduction-band spin-splitting in 
III-V QWs£ However, this is not the case for the L-point 
due to large energy spacing between the top valence band 
Ls' (X' , Y') and lower lying valence-band states. There- 
fore, one should again take into account the interband 
spin-orbit interaction. 

In the bulk, the spin-orbit coupling between the Lx 
and L31 bands is forbidden. However, the group C s has 
only two symmetry operations: identical and reflection 
in the (y', z')"P lane - The reflection yields (S'\U X >\X') = 
(S'\U y ' iZ '\Y') — 0, but the following matrix elements are 
nonzero: 



U{ = (S'\U y ,\X'), 
U' 2 = (S'\U Z ,\X'), 
U's = (S'\U X ,\Y>). 



(10) 



Introducing the matrix element of the bulk k • p- 
interaction V 1 = (U/mo) (S'\p x \X') , we obtain from 
Eqs. (JTJ and JTUI) for the i-states in a (001) QW 

2T" 

ATics = — ^rS{z-z i f)(U' 2 (j z ik x i+U' 1 ayik x i+U' 3 a xl (k y :)). 

(11) 



Here E' Q is the energy gap between the Lx and the L 3 > 
states and the angular brackets mean averaging over 
the quantum-confined state. As a result the component 
k y i = (k x + k y — 2fc z )/v / 6 reduces to (Ay) = (k x + k y )/^/6 
whereas the component k x i — (k x — k y )/\/2 remains un- 
changed. 

Inclusion of the spin splitting A of the valence band 
slightly modifies the coefficients in the above expression: 



AH Cs 



-S(z - z if )V 
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A/2 


E'o- 


hA/2^ 



UnO ' Z ' k x ' 



a x i(k 



v'l 



0~y> k x ' 



In (lll)-grown QWs the relevant symmetry is Cs v with 
the axis C3 and two more mirror reflection planes in ad- 
dition to the elements of C s group. The rotation yields 

U[ = -U's, U' 2 = 0, 

and we get from Eq. (|11|) the Rashba-like contribution 
for the (111) SiGe QWs with odd number of monolayers: 

2V'U[ 

AHc 3V = S ( z ' ~ z if) jp, {ox'ky, - o y ,k xl ). (12) 

If the number of monolayers is even (Dsd symmetry) then 
the inversion imposes the condition U[ = 0, and spin 
splitting is absent. 

Within the envelope function approximation, the elec- 
tron wave function tp satisfies the Schrodinger equation 
with the effective Hamiltonian 

H = H (k z ,k) + V(z) + AH. 

Here Ho is the bulk spin-independent Hamiltonian with 
k z = —id/dz, V(z) is the heteropotential, and the cor- 
rection AH is given by Eq. JHJ. 

Instead of solving the Schrodinger equation with the 
Hamiltonian H, one can equivalently find general solu- 
tions of the equations (Ho 4 V)tp — Eip within each ho- 
mogeneous layer and then apply the boundary conditions 



<p(z if -0) = VKf + 0), 



(13) 



2i . 



where the velocity operator v z — ftT 1 dHo/ dk z . 

In order to estimate the spin splitting one needs to go 
beyond the envelope function approximation. This can 
be done in the pseudo-potential or tight-binding model 
which yields the matrix Hi f in boundary conditions 11311 . 
The work on estimation of the interface-induced spin 
splitting in the microscopic tight-binding model is in 
progress. Here it suffices to note that the spin matrix 



4 



Hif in Eq. is non- vanishing if interatomic contribu- 
tions to -ffso either Hf. p are taken into account. 

In this paper, we have analyzed spin splitting of elec- 
tron states in SiGe heterostructures of all possible sym- 
metries. The absence of inversion center can be also 
probed by means of second-harmonic generation ex- 
periments ALifii An ideal Si m Ge n superlattice with odd 
n and m allows second-harmonic generation. It was 
demonstrated experimentally that Si m Ge„ super lattices 
with nominally odd and even n, m possess comparable 
second-harmonic conversion efficiency^ The weak non- 
linear response can be explained by antiphase micro- 
scopic domains shifted with respect to each other by one 
monoatomic layer along the growth direction [001]. 

Similarly, different domains are characterized by oppo- 
site signs of the linear in k spin-dependent matrix Hif 
in Eq. |J5J. The influence of this kind of imperfection on 
the D'yakonov-Perel' spin relaxation of the conduction 
electrons depends on the relation between the linear di- 
mension, Id, of a single domain and the electron mean 
free path, I, in the interface plane. In the case Id 3> I, 
the spin relaxation time of free carriers is the same as in 
a perfect structure: 

- = ^ (t p Tt (H? f )) [^{z lf )ip{z lf )} 2 . (14) 

Here t s is the relaxation time of spin z-component (for 
in-plane spin the corresponding time is 2t s ). In Eq. 1141) 



t p is the momentum scattering time, the angular brackets 
mean averaging over the carrier energy distribution, and 
ip(zif) is the interface value of the envelope function cal- 
culated at Hif = 0. In the opposite limiting case Id <C I, 
one has 



where t s is given by Eq. I|14f> . If the time t p is governed by 
scattering from the antiphase-domain walls (boundaries) 
one can use Eq. 1|14|1 for estimations of the spin relax- 
ation times. Thus, we conclude that, even if the overall 
symmetry of a SiGe heterostructure is due to the 
antiphase domains, the lack of inversion center within a 
particular domain leads to the spin dephasing according 
to Eqs. i|14[l and (|15|) . This opens a possibility to discuss 
the recent spin relaxation times measurements in these 
heterostructures 

In conclusion, we have shown that linear-in-fe spin 
splitting is present even in symmetrical SiGe QWs. It can 
be of Rashba, Dresselhaus, or both types. The splitting- 
is caused by anisotropy of chemical bonds at interfaces. 
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